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Abstract. The CAIRNS (Cluster And Infall Region Nearby Survey) project is a large spectro- 
scopic survey of the infall regions surrounding nine nearby rich clusters of galaxies. I describe 
the survey and use the kinematics of galaxies in the infall regions to estimate the cluster mass 
profiles. At small radii, these mass profiles are consistent with independent mass estimates from 
X-ray observations and Jeans analysis. I demonstrate the dependence of mass-to-light ratios on 
environment by combining these mass profiles with Two-Micron All-Sky Survey (2MASS) pho- 
tometry. Near-infrared light is more extended than mass in these clusters, suggesting that dense 
cluster cores are less efficient at forming galaxies and/or more efficient at disrupting them. At 
large radii, galaxy populations in cluster infall regions closely resemble those in the field. The 
mass-to-light ratio at these radii should therefore be a good probe of the global mass-to-light 
ratio. The mass-to-light ratio in the infall region yields a surprisingly low estimate of Q m ~ 0.1. 



1. Introduction 

The relative distribution of matter and light in the universe is one of the outstanding 
problems in astrophysics. Clusters of galaxies, the largest gravitationally relaxed objects 



in the universe, are important probes of the distribution of mass and light. Zwicky (1933) 
first computed the mass-to-light ratio of the Coma cluster using the virial theorem and 
found that dark matter dominates the cluster mass. Recent determinations using the 
virial theorem yield mass-to-light ratios of M/Lsj ~ 250/iAf Q /L Q (Girardi et al. 2000 
and references therein). Equating the mass-to- light ratio in clusters to the global value 
provides an estimate of the mass density of the universe; this estimate is subject to sig- 
nificant systematic error introduced by differences in galaxy populations between cluster 
cores and lower density regions fCarlbcrg et al. 1997[IGirardi et aT~ 2000l. Indeed, some 
numerical simulations suggest that cluster mass-to-light ratios exceed the universal value 
HDiaferio 19991 |Kravtsov fc Klypin 1999|lBahcall et al. 2000jl . 

Determining the global matter density from cluster mass-to-light ratios therefore re- 
quires knowledge of the dependence of mass-to- light ratios on environment. Bahcall et 
al. (1995) show that mass-to-light ratios increase with scale from galaxies to groups to 
clusters. Ellipticals have larger overall values of M/Lb than spirals, presumably a result 
of younger, bluer stellar populations in spirals. At the scale of cluster virial radii, mass-to- 
light ratios appear to reach a maximum value. Some estimates of the mass-to-light ratio 
on very large scales (>10/i _1 Mpc) are available (Bahca TTet al. 1995|> . but the systematic 
uncertainties are large. 

There are few estimates of mass-to-light ratios on scales between cluster virial radii 
and scales of 10/i _1 Mpc (Rines et al. 2000, Rines et al. 2001a, Biviano & Girardi 2003, 
Katgert et al. 2003, Kneib et al. 2003, Rines et al. 2004, Tully 2004 and references 
therein). On these scales, many galaxies near clusters are bound to the cluster but not yet 
in equilibrium IjOunn & Gott 1972|l . These cluster infall regions have received relatively 
little scrutiny because they are mildly nonlinear, making their properties very difficult to 



1 



2 



K. Rines 



predict analytically. However, these scales are exactly the ones in which galaxy properties 
change dramatically (e.g., Ellingson et al. 2001, Lewis et al. 2002, Gomez et al. 2003, Treu 
et al. 2003, Balogh et al. 2004, Gray et al. 2004). Variations in the mass-to-light ratio with 
environment could have important physical implications; they could be produced either 
by a varying dark matter fraction or by variations in the efficiency of star formation 
with environment. In blue light, however, higher star formation rates in field galaxies 
could produce lower mass-to-light ratios outside cluster cores resulting only from the 
different contributions of young and old stars to the total luminosity (Bahca TTet al. 2 000 
|Tully 2004| ). 

Galaxies in cluster infall regions produce sharp features in redshift surveys. Early in- 
vestigations of this infall pattern focused on its use as a direct indicator of the global 
matter density Q m . Unfortunately, random motions caused by galaxy-galaxy interactions 
and substructure within the infall region smear out this cosmological signal (Diaferio & 
Geller 1997, Vedel & Hartwick 1998). Instead of sharp peaks in redshift space, infall 
regions around real clusters typically display a well-defined envelope in redshift space 
which is significantly denser than the surrounding environment (IRines et al. 20031 here- 
after Paper I, and references therein). 

Diaferio & Geller (1997) and Diaferio (1999) analyzed the dynamics of infall regions 
with numerical simulations and found that in the outskirts of clusters, random motions 
due to substructure and non-radial motions make a substantial contribution to the am- 
plitude of the caustics which delineate the infall regions. Diaferio & Geller (1997) showed 
that the amplitude of the caustics is a measure of the escape velocity from the cluster; 
identification of the caustics therefore allows a determination of the mass profile of the 
cluster on scales < 10/i _1 Mpc. 

Diaferio & Geller (1997) and Diaferio (1999) show that nonparametric measurements 
of caustics yield cluster mass profiles accurate to ^50% on scales of up to 10 hT x 
Mpc. This method assumes only that galaxies trace the velocity field. Indeed, simu- 
lations suggest that little or no velocity bias exists on linear and mildly non-linear scales 
(Kauffman n et al. 1999al|Kauffmann et al. 1999b|l . The caustic method has been applied 
to systems as large as the Shapley Supercluster jRcisenc gger et al. 2000| ) and as small as 
the Fornax cluster fDrinkwa ter et al. 2001(1 as well as to many nearby clusters (Paper I). 
Biviano & Girardi (2003) applied the caustic technique to an ensemble cluster created 
by stacking redshifts around 43 clusters from the 2dF Galaxy Redshift Survey. Rines et 
al. (2000) found an enclosed mass-to-light ratio of M/Lr ~ 300h within 4 h~ 1 Mpc of 
A576. Rines et al. (2001) used 2MASS photometry and the mass profile from Geller et 
al. (1999) to compute the mass-to-light profile of Coma in the K-band. They found a 
roughly flat profile with a possible decrease in M/Lk with radius by no more than a fac- 
tor of 3. Biviano & Girardi (2003) find a decreasing ratio of mass density to total galaxy 
number density. For early-type galaxies only, the number density profile is consistent 
with a constant mass-to-light (actually mass-to-number) ratio. 

Here, we calculate the infrared mass-to-light profile within the turnaround radius for 
the CAIRNS clusters (Paper I), a sample of nine nearby rich, X-ray luminous clusters. We 
use photometry from 2MASS, the Two Micron All Sky Survey ( Skrutski e et al. 1997i and 
add several new redshifts to obtain complete or nearly complete surveys of galaxies up 
to 1-2 magnitudes fainter than (as determined by Cole et al. 2001 and Kochanek et 
al. 2001). Infrared light is a better tracer of stellar mass than optical light; it is relatively 
insensitive to dust extinction and recent star formation. Despite these advantages, there 
are very few measurements of infrared mass-to-light ratios in clusters ( Tust in et al. 200 ll 
IRines et al. 2001^1 ILin et al. "20031 IRines et al. 200411. 
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Figure 1. Redshift versus projected radius for the combined CAIRNS cluster. The solid lines 
indicate the caustics and the errorbars show 1-a uncertainties. 

2. Observations 

The CAIRNS clusters are 8 of the 14 Abell clusters that are: nearby (cz < 15, 000 km s^ 1 ), 
Abell richness class R^l, X-ray luminous [Lx > 2.5 x 10 43 /i _2 erg s _1 ) galaxy clusters 
with declination 6 > —15°. Between 1997 and 2003, we collected 5607 redshifts (both 
new and remeasured) in large areas around nearby clusters with the FAST spectrograph 
on the FLWO 1.5-m telescope in Arizona. We targeted galaxies within a projected ra- 
dius of ~10ft, Mpc of the clusters, selecting targets first from POSSII 103aE plates and 
later from 2MASS when it became available. The redshift catalogs in Rines et al. (2004) 
complete K s band magnitude- limited samples extending 1-2 magnitudes fainter than the 
characteristic magnitude MS- . These samples include ^60-85% of the total light in the 
clusters and their infall regions. 

3. Mass Profiles 

Figure ^ shows the redshifts of galaxies surrounding nearby clusters as a function of 
projected radius (normalized to r 2 oo, the radius within which the average density is 200 
times the critical density). The caustic pattern (a dense envelope in phase space with 
well-defined edges) is evident in the combined cluster as well as in each of the 9 clusters. 
Using the phase space distribution of galaxies in cluster infall regions, we apply the 
kinematic mass estimator of D99 to these clusters. The resulting mass profiles agree well 
with NFW (Navarro et al. 1997) and Hernquist (1990) models, but exclude an isothermal 
sphere. 

These mass profiles agree with X-ray masses at small radius as well as virial masses 
(after correction for the surface pressure term) at slightly larger radii plines et al. 20021 
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Figure 2. if s -band luminosity functions for the CAIRNS clusters. The left panel shows galax- 
ies in the virial regions (projected within r2oo); the right panel shows galaxies in the infall 
regions (projected outside r2oo)- Solid lines show the shape of the 2dF/2MASS its-band LF 
with arbitrary normalization. 



Paper I) . The latter is primarily a consistency check as the caustic technique utilizes the 
same kinematic data as the virial theorem. This consistency can be further demonstrated 
with the velocity dispersion profiles (see Paper I and Rines et al. 2004). 



4. Near-Infrared Luminosity Functions 

When using the mass-to-light ratio in clusters to estimate il m -, one must determine 
whether the luminosity functions of field and cluster galaxies differ significantly. Given 
the well-known morphology-density relation, it is possible that the two LFs differ signif- 
icantly. We use 2MASS photometry to determine the near-infrared luminosity functions 
for the CAIRNS clusters and infall regions ( Rines et al. 2004(l . The cluster and infall re- 
gion LFs are very similar to each other (Figure |2J) and to the field LF for relatively bright 
galaxies (Mpc s < MJ- +2). Because of this similarity we can correct for the luminosity in 
faint galaxies using the field LF. Our redshift surveys include 60-85% of the total light. 



5. Near-Infrared Mass-to-Light Profiles 

Both the surface number density and surface luminosity density profiles of cluster /infall 
region members are more extended than the mass profiles of Paper I. Figure shows the 
mass-to-light profiles of the CAIRNS clusters. The mass-to-light profiles are either flat 
or show a decreasing M/Lr with increasing radius. The mean mass-to-light ratio within 
r2oo is a factor of 1.8 ± 0.3 larger than the mean value outside r2oo- 

The decreasing mass-to-light profiles could be caused by gradients in stellar populations 
with radius; such effects have been invoked to account for a similar result for B band 
mass-to-light profiles in simulations (Bahc aTTet al. 2 000). In K band, however, changes 
in stellar populations with radius are expected to change the mean stellar mass-to-light 
ratio by < 20%. Thus, the CAIRNS mass-to-light profiles provide tentative evidence for 
variations in the efficiency of galaxy formation and/or disruption. Environments with 
higher virial temperatures (like cluster centers) are more efficient at disrupting galaxies 
and/or less efficient at forming them. 
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Figure 3. Mass-to-light ratio in K s band as a function of radius for the CAIRNS clusters. 
The solid lines show the caustic mass profile M(< r) divided by the projected luminosity pro- 
file Lk s {< R P ) and the shaded regions indicate the associated 1-cr uncertainties. The open 
squares show the mass-to-light ratio M(r,r + dr)/ Lk s (R p , R p + dR p ) in radial shells. The 
dash-dotted line shows the projected best-fit Hernquist mass profile Mh{< R p ) divided by 
Lk s (< R p )- The solid stars show the mass-to-light ratio calculated using the X-ray tem- 
perature and the mass-temperature relation to estimate the mass. The lower point shows 
M x (< r 500 )/L Ks (< R500) and the upper point shows M x (< Rsoo)/Lk s (< R500) assuming 
Mx(< -R500) = 1.3Mx(< )~5oo) as is true for an NFW mass profile with c~5. 



A related trend has been noted by Lin et al. (2003, 2004), who show that K-band mass- 
to-light ratios within r§oa increase with increasing cluster mass. This result indicates 
that more massive clusters (with larger virial temperatures) have less efficient galaxy 
formation and/or more efficient galaxy disruption. We confirm this trend in Rines et 
al. (2004). Because cluster infall regions should be composed of galaxies inhabiting less 
massive systems | Rines et al.(2001b| IRines et al. 2002|l and/or regions with lower virial 
temperatures, the above trend predicts that cluster infall regions should have smaller 
mass-to-light ratios than virial regions, consistent with the CAIRNS results. 

Figure 0] shows the mass-to-light profile of the combined CAIRNS cluster (normalized 
to unity at r2oo)- The mass-to-light profile clearly decreases with radius. The combined 
cluster should be less susceptible to substructure than the individual clusters. Figure 0] 
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Figure 4. Mass-to-light ratio in K s band as a function of radius for the combines CAIRNS 
cluster. The thick and thin solid lines show the caustic mass profile M(< r) divided by the 
projected luminosity profile Lk s {< Rp) and the associated l-cr uncertainties. The open squares 
show the mass-to-light ratio M(r, r + dr) / Lk s {Rp, Rp + dR p ) in radial shells. 

shows that the mass-to- light ratio in shells decreases by about a factor of 2 from the virial 
region to the infall region, consistent with the results for individual clusters. Taking the 
mass-to-light ratio in cluster infall regions as an estimate of the global value and the 
SDSS luminosity density extrapolated to K s band IjBlanton et al. 2003(1 . we estimate 
f2 m = 0.10 ± 0.02 (statistical). This estimate is somewhat lower if we use the K s band 
luminosity density of either Cole et al. (2001) or Kochanek et al. (2001). We discuss 
potential systematic effects in Rines et al. (2004). Most of these effects would flatten the 
observed mass-to-light profiles relative to the true profiles, suggesting that the observed 
decrease in mass-to-light ratio with radius is real. 

6. Conclusions 

Cluster infall regions contain more galaxies than their virial regions. If currently pop- 
ular cosmological models are correct, the mass in infall regions will eventually accrete 
onto the parent clusters, and their final masses will increase by a factor of about 2. 

Near-infrared luminosity functions depend only weakly on environment, at least at 
the bright end. We show that near-infrared light is more extended than mass in cluster 
infall regions, suggesting environmental dependence of the efficiency of galaxy formation 
and/or disruption. If more efficient galaxy disruption is responsible, intracluster stars 
might be a significant component of stellar mass in clusters. The mass-to-light ratios in 
infall regions suggest a low fi m ~ 0.1. Future work is needed to determine the significance 
of the conflict of this result with the currently favored Q, m ~ 0.3. 
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